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We give a universal recipe for constructing nonlinear entanglement witnesses able to detect non- 
classical correlations in arbitrary systems of distinguishable and/or identical particles for an arbitrary 
number of constituents. The constructed witnesses are expressed in terms of expectation values of 
observables. As such they are, at least in principle, measurable in experiments. 
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■ I. INTRODUCTION 

X) ■ 

I ' Nonclassical correlations among subsystems of a composite quantum system, known as entanglement, can be easily 
, characterized mathematically but, at least in the case of mixed states, in a rather ineffective way. In general, given 
' a mixed state of a composite system it is hard to decide whether it is separable (nonentangled) or not with respect 
\ to a given partition of the whole system into subsystems. In the case of many (more than two) subsystems even 
I— I ■ determination of separability of a pure state might pose a computational problem. 

' It is thus desirable to construct a 'measure of entanglement', i.e. a function from the set of quantum states into 
real numbers which vanishes on separable states and takes nonzero, say positive, values for the nonseparable ones. If 
we impose some further, natural and reasonable conditions to be fulfilled by such a function, we may obtain a useful, 
quantitative characterization of 'the amount of entanglement' in a given state. Evidently such a quantitative measure 
^ is not unique [l| and the choice of a particular measure is dictated by a particular case we want to analyze. In the 
following we will mostly focus our attention on discriminating between separable and nonseparable states. In terms 

■ of an entanglement measure it means that we are interested only if it takes zero or non-zero value on a state under 
^ ■ investigation. 

^ . To be of practical use in experiments, an entanglement measure should be measurable, i.e. it should be possible 
' to design an experiment in which a value of the measure for a particular state can be established. Since in quantum 
^-H , mechanics we can measure only observables, a measurable measure should be given as an expectation value of a 
C\| Hermitian operator calculated in the investigated state. It is, however, straightforward to check that there is no 
] observable such that its expectation values vanish exactly on pure separable states |2|]. Such known measures of 
' entanglement like, e.g., Wootter's concurrence can be expressed in terms of expectation values of antilinear operators 
' and as such are not directly measurable. It is however possible to find bilinear (i.e. acting on two copies of a 

. state) Hermitian operators for which the condition of vanishing expectation values are fulfilled only by pure separable 
states 

In the following we show how to construct such measures for pure states in the two-partite as well as many-partite 
systems in an 'algorithmic' way. Moreover the constructed measures allow to estimate from below the amount of 
. , entanglement for mixed states and thus provide some effective (although not always decisive) criteria of entanglement. 
• One of our main points is to stress the unifying character of the presented approach, allowing for a uniform treatment 
of arbitrary number of distinguishable and identical particles (bosons and fermions) . 
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II. PURE SEPARABLE STATES 



As it is customary in quantum information theory, for which entanglement is one of the most important resources, 
we will investigate quantum systems in finite-dimensional Hilbert spaces (one should think about various spin or spin- 
like system, multilevel atoms etc.). Thus with a quantum system we associate a Hilbert space H isomorphic to C^. 
Customarily, vectors from H are called (pure) states of the system. One should, however, keep in mind that in order 
to give a probabilistic interpretation to amplitudes we normalize states to unit norm (length) and moreover disregard 
an irrelevant phase of the vector. It is thus more appropriate to think of states as points in the projective space 
P('H). Alternatively and equivalently, it is convenient to treat a pure state as a one-dimensional projection operator 
:= freeing ourselves from the normalization and phase problems and unifying the treatment of pure 

and mixed states by identifying both with positive-definite operators on H with unit trace, where pure states are 
distinguished hy — p exhibiting their projective character. 
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For composed systems the Hilbert space Hcomp is a tensor product of the Hilbert spaces of subsystems, 'Hcomp = 
^1 (8) 'H2 (8) • • • <8) 'Hn- The states represented by shnple tensors, i.e. 'Hcomp 9 |$) — \4'i)® 102) |0n) are called pure 
separable states. Identifying a pure state with a one-dimensional projection operator we call such a projection separable 
if it projects on the direction of a separable vector. We then define separable mixed states as the convex hull of the 
separable projections. It is easily seen that this definition of separable states coincides with the commonly adopted one 
- a state is separable if it is a combination of product states with positive coefficients, p — J2iPiPi ^ pf ■ pf, Pi - 
a state on Hi, Pi > 0. The states that are not separable are called entangled. From the point of quantum information 
technology the really interesting states, useful, e.g., in secure and more effective transmission of information, are 
entangled states, and that is why effective discrimination between them and separable ones is of importance (see, e.g., 

H). 

In the case of bipartite systems an arbitrary pure state can be represented by a vector 



i,3 



(1) 



where {\ei)}!f^i and {\fj)}jLi are some orthonormal bases in Hi (of dimension N) and (of dimension M), and 
c — (cki) is a complex matrix. By unitary changes of the bases jci) = J2k Uki\e'k), \fj) = J2i which amounts 

to transformation c 1— > UcV'^ '■— c' , with ^ standing for transposition, one can bring ([T]) to its Schmidt form. 
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Afcl4)®l/fc), Afc>0. 



(2) 



A pure state is nonentangled if and only if only one of its Schmidt coefficients Afc does not vanish. The coefficients 
Afe are positive real numbers squares of which are non-zero eigenvalues of c^c (or, equivalently cc^), hence are easily 
calculable for a given | ■(/'), 

The definition of separability in terms of simple tensors lacks sense for systems of identical particles when states 
must be symmetric or antisymmetric with respect to relabeling of subsystems (particles). Indeed due to the 
(anti)symmetrization nearly all state vectors do not have the product form; the only exception is the state of bosons 
each occupying the same single particle state). For example, the simplest state vector of two identical fermions reads 
as 1$) — |0i) (g) 102) — 102) <8) |0i) |0i) A 102) which according to the usual definition is entangled (33| . 

In a way, states of identical particles exhibit some a priori correlations and only additional amount of correlation 
should be classified as entanglement. Several, not necessarily equivalent ways were proposed to identify and quantify 
this phenomenon. In a correlation measure for states of two undistinguishable fermions was proposed. It 

is constructed in analogy with the distinguishable particles case by employing algebraic properties of the coefficient 
matrix in the expansion of a state in terms of basis states. To this end one observes that a pure states of two 
undistinguishable fermions in an n-dimensional single particle space H can be written as 



(3) 



where w — (wij) is a complex antisymmetric matrix fulfilling the normalization condition ti{w''w) = ^. Here are 
fermionic creation operators and |0) is the vacuum state. A unitary transformation U in the single particle space 
changes w to w' = UwU'^ . By choosing an appropriate unitary U we can transform w to its canonical, block-diagonal 
form: 



= diag[Zi,...,2'r,Zo], Z.i = 











> 0, 



(4) 



where 2r is the rank of w and Zq is the null matrix of dimension {n— 2r) x (n~2r) [13. The squares of the coefficients 
Zi are eigenvalues ofww'' , hence again are easily calculable for a given w. A state is, by definition, nonentangled if only 
one of Zi does not vanish, i.e. only a single, elementary 2x2 Slater determinant appears in the canonical decomposition 
(|4]). In other words a state is nonentangled if it can be written as {ip 



— /i /2 I*-*) where |0), i = 1, . . . , n form an 
orthonormal basis in the single particle Hilbert space, i.e. \^) is the antisymmetrization of a product state. Measuring 
quantum correlations by the number of nonvanishing terms when expanding a wave function in terms of elementary 
Slater determinants was proposed earlier [isl , [l9| , employing various quantitative characterizations of the number of 
non-zero terms. 

The above idea was thoroughly investigated in [l^,[l^[23|. Generalizing the two-particle case, the relevant definition 
of pure nonentangled state of n indistinguishable fermions (we will use the notion of 'nonentanglement' rather than 
'separability' which in this context lacks its semantic sense), can be shortly summarized as follows. A fermionic state 



3 



is nonentangled if and only if it can be written as f} f J • • • /^|0), i.e. it is the complete antisymmetrization of a product 
state. This coincides with definitions proposed in |2ll [2^ , where conclusions were reached by departing from slightly 
different starting points. 

Similar ideas can be applied to bosons [13, Hll, [2^ . A general, two-particle state in an n-dimensional single particle 
space, 

k)= E%-^M|0), (5) 

where &J are boson creation operators and v = (vij) is a complex symmetric matrix, can be transformed to a 
canonical form with Vij = ZiSij by a unitary transformation in the single-particle space, amounting to v U vU'^ 
on the level of the coefficient matrix (33 | . Consequently a pure n-boson state is nonentangled if it can be written as 
\v) = b\ ■ ■ ■ b\\0) . It should be pointed that in [2lj and (see also[25|), in contrast to [10, 113|, a slightly broader 
definition of pure nonentangled boson states was proposed. In addition to the identified above, as nonentangled are 
also treated states which in some basis can be written a.s b\ ■ ■ ■ b\f\0) , where all states bk\0) are orthogonal. We would 
like to make two comments concerning this point. Firstly, the extended class of states can be also easily described 
using methods proposed in this paper (as orbits of unitary groups - see below), although admittedly, characterization 
in terms of vanishing expectation value of some bilinear operator needs more efforts which we postpone to other 
occasion. Secondly, the real meaning of entanglement becomes important in particular experiments. Whether it can 
be exhibited (especially for indistinguishab le p articles) depends strongly of what and how we measure states in an 
experiment [Ig, thus, e.g., the authors of [13| mention that states identified by Li et al. [ll| as nonentangled can 
be used for quantum teleportation, and as such should be treated as quatally correlated (entangled). In the present 
paper we would like to retain the more restrictive definition of bosonic nonentanglement, leaving, as already stated, 
the more relaxed one for the topic of future investigations. 

We would like to conclude this section by stressing that there are no analogs of the Schmidt and Takagi decompo- 
sitions in multipartite cases, therefore establishing nonentanglement demands considerably more elaborate methods 

MM- 

III. ACTIONS OF UNITARY GROUPS IN COMPOSITE SYSTEMS SPACES 

In order to achieve the goals of characterizing separability via observables let us look at the problem from a 
slightly different point of view. For the moment we restrict the attention to bipartite systems. It is obvious that 
separability of a pure state of distinguishable particles does not change if we individually evolve each subsystem in 
a quantum-mechanically allowed way i.e. via a unitary transformation, |02) Ui\(j)i) U2\4>2)- In the case 
of indistinguishable particles the same is true if we perform the same unitary transformation U in each one-particle 
space, l^i) A |02) A ?7|02) or j^i) V |(/>2) U\(f>i) V U\(j)2) (where V denotes the symmetrized tensor product). 

In fact we used this invariance to transform pure states to their canonical Schmidt or Takagi forms in the preceding 
section. We may assume that the matrices have determinant one, so they are elements of special unitary groups, since 
the phase of the state does not play any role. The action of the special unitary group on vectors from a Hilbert space 
T-l translates in a natural way to an action on the projective space P{H). If we denote by [v] the point in F{H) (the 
direction of \v), or in other words the complex line through \v)) we have, by definition, U[v] = [U\v)]. This action of a 
unitary group on P(?^) is transitive, i.e. any two points in P(?^) can be connected by some unitary transformation. A 
straightforward conclusion is that any two nonentangled states in spaces of states of bipartite systems are connected by 
the above described unitary action of the direct product of two unitary groups in the case of distinguishable particles 
and a single unitary group for fermions and bosons. In more technical terms nonentangled state in all cases form a 
single orbit of an action of the group SU{N) x SU{M) in P('Hi ® H2) for distinguishable particles or SU{N) in the 
case of bosons or fermions in P{H V H) or P{H A H), respectively. 

In order to identify uniquely the orbit in question we have to make a short excursion in the theory of representation 
of semisimple Lie groups (the special unitary group SU (N) is simple and the direct product of simple groups, such as 
SU{Ni) X SU{N2) X • • • X SU{Nk) is semisimple 27]). Let if be a semisimple real group and G its complexification, 
and denote by g the Lie algebra of G. For K = SU{N) we have G — SL{N,C) (special complex linear group in N 
dimensions) and g = s[jv(C) (the algebra of complex N x N matrices with vanishing trace). As a basis in sIn we 
can choose N — 1 independent traceless diagonal matrices Hk which span its maximal commutative subalgebra and 
the matrices Xij having a single nonvanishing element on in the position. The commutation relations among 
Hk and Xij read as: [Hk^Xij] — aij{Hk)Xij where is some linear function on the set of diagonal matrices. An 
analogous construction exists for an arbitrary semisimple complex Lie algebra. We can distinguish in it a maximal 
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commutative subalgebra t of dimension r (called the rank of the algebra) and one-dimensional subspaces (root spaces) 
Qa (spanned by X^) such that 

[H,X]^a{H)X, Hei, Xegc. (6) 

If we choose a basis {Hk}l^i in i, the algebra is uniquely determined by the set of vectors a — {a{Hi), . . . ,a{Hr)) 
where a runs over all different root spaces. The linear form a is called a root of g, and the element Xa corresponding 
to the root a is a root vector. There always exists a natural symmetry: to each a corresponds —a, as in the described 
above case of slAr(C) treated as the algebra of N x N traceless matrices, to each Xij with, say, i < j, i.e. with a 
single nonvanishing element in the upper-right triangle, there corresponds Xji living in the lower-left triangle of the 
matrix for which the commutators with Hk have the opposite sign to those of the commutators of Xij. It means 
that to characterize an algebra we need only half of the root vectors - the 'positive' ones corresponding e.g., to upper 
triangular matrices which we will call positive root vectors and the corresponding roots a - the positive roots (a > 0). 

The groups K and G, as well as the Lie algebra g can be represented irreducibly in spaces of different dimensions, 
i.e. to each element of the group or algebra there corresponds a linear operator acting in some Hilbert space, say, C^, 
such that the group multiplication and the Lie bracket (commutator) in the algebra are preserved. Hence, denoting 
by Tr{X) the representative of the Lie algebra element X we have from ([6]) 

[Tr{H),niX)]^aiH)7r{X), H e t, X G (7) 

Since the operators {Hk}k^i commute, the same is true for their representatives {TT{Hk)}%^i- It follows that Tr{Hk) 
have common eigenvectors. For each irreducible representation there exists a unique (up to a multiplicative constant) 
eigenvector \vmax) of all T:{Hk) which is annihilated by all the representatives of the positive roots 

TT{Hk)\Vynax) = '>^k\Vmax) , fc = l,...,r, ■n{Xa)\Vmax) , « > 0. (8) 

The vector \vmax) is called the maximal weight- vector. An irreducible representation of g (and, in consequence, of G 
and K) is uniquely determined by the eigenvalues Afe which we cast in a r-component vector A = (Ai, . . . , A,.). 

It is now easy to identify the orbit of SU{N) x SU{M) in P(Hi ® H2) of the nonentangled states as the orbit 
through the maximal weight vector for the representation of SU{N) x SU (M) in Hi (81 H2 — (X) C*^. Indeed, since 
the action of SU{N) on P(C^) is transitive, each point [v] is on the orbit {[J7|vmaa;)] : U £ SU{N)}. Thus each 
nonentangled state is on the orbit of SU{N) x SU{M) of the state [uimaa;] := [vmax Umax] where \vmax) and \u„iax) 
are the highest-weight vectors of the representations of SU{N) in and SU{M) in C-'^ . On the other hand, it is 
easy tho see that \wmax) is exactly the highest-weight vector of the irreducible representation SU{N) x SU{M) in 
Hi (E) H2 — ® C*^. Similar considerations show that analogous statements are true in the cases of many particles 
with and/or without assumptions about their distinguishability. 

Summarizing, the nonentangled states form the orbit in the projective space P('Hcomp) through the highest-weight 
vector of an irreducible representation of a unitary group K in Hcomp- For ^-distinguishable particles K — SU{Ni) x 
■ • • X SU{Nn) and Hcomp = "Wi ^E) ■ ■ ■ <E) Hn, dimHk = N^, whereas for indistinguishable particles K = SU{N) and 
'^comp = HV ■ ■ ■ V H (for bosons) or Hcomp — H /\ ■ ■ ■ AH (for fermions), diuiH — N. 

There exists a nice and simple method of characterizing the orbit through highest weight vector of an irreducible 
representation of a semisimplc group (28j . First we define a second-order operator 

r 

C2 := {AXc.)7r{X^o.) + 7r(X_„)^(X,)) + ^ ^{Hi)^{H,). (9) 

Q>0 i=l 

It is called the (second-order) Casimir operator; using ([7]) one shows that it commutes with all operators of the 
representation of the algebra, and thus (for an irreducible representation) is proportional to the identity operator I 
|35| . In fact one can prove that C2 = (A, A -I- 25)1 , where 5 = ^ X]a>o half-sum of the positive roots and ( , ) 

is the Euclidean scalar product in the space of r-dimensional vectors to which A, a, and 5 belong p9| . 

A vector belongs to the highest- weight orbit of the irreducible representation tt with the highest weight A of a 
semisimple group G if and only if 

L|^)(X|V) = (2A + 2<5,2A)|V')0|^), (10) 

where 

r 

i = C2«)/ + /®C2 + 2^ (7r(A:„)®7r(A:_„) + 7r(X_„)(g)7r(X„)) + 2 ^ 7r(iJ,) «) 7r(i?0- (H) 

Q>0 1=1 

Relevant properties of the operator L for the cases of two distinguishable particles, fermions, and bosons are 
calculated in the Appendix. For general cases of n particles, distinguishable or not, similar calculations can be also 
performed (with considerably more effort). 
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IV. MEASURABLE ENTANGLEMENT MEASURES 



The operator L is Hermitian and positive semidefinite and its largest eigenvalue equals Imax = (2A + 2^,2A). 
Hence A := ImaxI — L \s positive semidefinite and its expectation value ® (V'l^lV') ® IV') vanishes exactly for a 
nonentangled This is the desired characterization of the nonentangled states in terms of an expectation value of 
an observable, i.e. A can be regarded as a 'nonlinear' entanglement witness. 

The operator A acts in the tensor product Jicomp T^comp- There is a natural isomorphism (the Jamiolkowski 
isomorphism [30l |) between the set of operators on T-icomp 'Hcomp (such as the operator A) and the set of operators 
acting on density matrices (or generally linear operators) on T-Lcomp given by 

A(p) = tri((p^®/)A), (12) 

where denotes the transposition of a matrix and tri is the trace over the first factor of the tensor product 'Hcomp ® 
"Hcomp- The most important feature of p2)) is, for our purposes, that for positively semidefinite A the operator A is 
completely positive. The actual definition of complete positiveness is not important here - what is crucial is that it is 
equivalent to the fact that the action of A can be expressed in the so-called Kraus form [l| 

s 

A(p)-^T^pr;. (13) 
The operators can be expressed in terms of eigenvectors of A. Let 

3 

^ = 5]i^^|«M)(i'Ml (14) 

be the spectral decomposition of A and s the rank of A. Since A is positive semidefinite we have t'^ > and defining 
\Wfj.) = we obtain 

s 

^ = El^/^)KI- (15) 

Let us now choose in Hcomp an orthonormal basis je^), i = I, . . . ,N = dim Hcomp, and expand the imnormalized 
eigenvectors \w^) of A in this basis, 

AT 

K)= E (16) 

It is now easy to find that 

N N N 

Another straightforward calculation reveals the matrix elements of A in terms of (sij 

s 

{i;2^HA\i^i^i^3)^Y.(i'2\Tf.m){ri\Tlm- (18) 

In particular, with the help of the expectation value of A we can construct a pure-state entanglement measure 
(generalized concurrence) 



1 /2 

{,p):={^^,p\A\^^,py/^ ^ (^Y.\{^\T^\r)\') , (19) 



CA 



with the property CAicttp) — |c«Pca(V') [1, 
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Since the Kraus operators are calculable from the eigenvectors of A and the latter are the same for L which 
differs from A by a multiple of the identity operator, we can find them by decomposing the image of L under the 
Jamiolkowki isomorphism into the Kraus form. When calculating from eigenvectors of L we should disregard 
those which correspond to the maximal eigenvalue Imax of L since they belong to the zero eigenvalue of A and as 
such do not contribute to the Kraus decomposition of its image under the Jamiolkowski isomorphism [cf. (|14ll5p ]. In 
fact, we may disregard even more of the Kraus operators calculated from L. Indeed due to the appearance of the 
complex conjugate vector lip*) in (jl9p . the terms {iplT^lip*} vanish for antisymmetric T^, hence for our purposes only 
symmetric are of interest when determining ca- For example, as shown in the Appendix, only the subspace 
of Hcomp 'Hcomp = 'Hi ^ 7^2 'Hi (8) "^2 obtained by antisymmetrizing separately the "Hi factors and the 7^2 ones 
produces relevant Kraus operators and A can be thus chosen as the projection on Tig fSG']. Similar analysis may be 
also performed in the case of bosons and fermions. In both cases A is a projector on subspaces and described 
in details in the Appendix. 



V. MIXED STATES 



Every mixed state p can be decomposed as a convex combination of pure states 

K 



P 



k=l 



which, with \(j)k) — v^Pfc|'0)fej can be rewritten as 



K 



p = (21) 



fc=i 

A particular example of ([20)) is provided by the spectral decomposition of p, 

R 



P = ^rk\'nk){'nkl p\Vk) ^ rklvk), R^rankp, (22) 
fc=i 



or, equivalently, as in (|2T|) 

R 



P = J2\^''){U \ik)^V^\rik)- (23) 
fc=i 



In fact all other decompositions ([2T|) can be obtained from (|23| via 

R 



m=Y.^kAQ, (24) 



where V is a. K x R matrix fulfiUing V^V = I "s^. 
We may now define for a particular decomposition ((2T 



K 

CA({0fc}) (25) 
fc=i 

where ca{iP) is given by ([T9|. It is now obvious that if minimizing ca({&}) over all decompositions ([21]) of p gives 
zero then p is nonentangled since it has a decomposition into a convex combination of nonentangled pure states. We 
obtain in this way a well defined measure of entanglement [3| for mixed states 

ca(/o) = min^Ca(0fe), (26) 
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where the minimum is taken over all decompositions ((20|) . Using ([T8| and ([24)1 we obtain further 
Ca{p) = min^CA(0fc) = min^(0fe (8)(/)fc|A|0fe ® (/)fc)i/2 

fc k 

= min^ (yfe*14* 14;Vfe„(e, $5 ^,\A\^i ® 
fc 

= min^ iv^^V^^VkiVkm Y.^mCj){Ci\Tl\^n ' ' 



fc At k fi 



fcfc 



1/2 



where are r x r matrices, 



(27) 



(28) 



constructed from easily calculated ingredients: the eigenvectors of p and the Kraus operators of A (or, equivalently, 
L), and the minimum is taken over all V fulfilling V^V. 

We may further reduce the complexity of minimization [BJ by employing the Cauchy-Schwartz inequality 



With — \(y*T^V^)kk\ and J^fj. VmU — 1. We obtain thus 



(29) 



CA 



{p) > min5]5^2/^|(T/V^yt)^^| > ^^J- 



k fJ. 



kk 



(30) 



where we used \z^ \ > 2^|. The minimization over V can be now performed explicitly giving, 



ca(p) > max ^ 0, Ai — ^ Aj 



(31) 



where A^ are the singular values of the matrix T = VfiTfj.- The matrix T still depends on the parameters which 
can be chosen in an arbitrary way under the condition = 1, leaving a large freedom to construct lower bounds 

for ca(p). 



VI. SUMMARY AND CONCLUSIONS 



We have presented a method of discriminating pure nonentangled states for multipartite systems. The method is 
universal - it applies, at least in principle, to systems with arbitrary number of distinguishable as well as undistin- 
guishable particles. Let us point at some other advantages of the proposed approach 

• The defined measure of entanglement is expressed in terms of a Hermitian, albeit bilinear operator - "a nonlinear 
entanglement witness" , and as such is, in principle, a physically measurable quantity. 

• The method, based solely on representation theory, can be easily adapted to more complicated situations, e.g. 
systems consisting of mixtures of bosons and fermions. 

• Calculation of the generalized concurrence is made in an algorithmic way consisting of few steps: 1) identification 
of a relevant group of local transformation and its representation, 2) calculation of the Lichtenstein's operator 
L given in terms of the operators of the Lie algebra of the group, 3) identification of relevant Kraus operators 
of the image of L under the Jamiolkowski isomorphism. To perform the last step one looks for the symmetric 
Kraus operators which are obtained from eigenvectors of L not belonging to its largest eigenvalue. The latter 
can be explicitly calculated from data about the group and representation in question. 

The pure state generalized concurrence constructed here can be used as a basis for effective estimates of mixed-state 
entanglement. One example of such estimate has been presented in the last section of the paper. 
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VIII. APPENDIX 

Let us choose the following basis in sl]\[{C) 

X^J = mi i,j = l,...,N, (32) 

Hi ^ ^J_-^ (^^|fc)(fc|-/|/ + !)(/ + 1|^ =^a;fc|fc)(fc|, l^l,...,N-l. (33) 

The normalization of the basis elements was chosen to have trHf = 1 = tr XijXji. The positive roots correspond to 
i <j in (|32|) . 

Short calculations show that 

tlt\ (34) 

Let us start with distinguishable particles. For simplicity we assume that Hi = H2 ='■ H, dimH = N. Remember 
that despite this we consider distinguishable particles, so we represent SU (N) x SU{N) on Hcomp — 'Hi®'H2 — 'H®'H. 
The corresponding representation of the Lie algebra s[Ar(C) ® slAr(C) is given by 



7r((A,0)) = 7r((0,B)) (35) 



hence 



7V-1 



Q>0 j=l 

Using the explicit representation ([5 ^ - p3p we obtain 

C2=[l-^y^I. (37) 

Let \ijkl) := G 'Hcomp®'Hcomp = 'H®'H®'H®'H and let us introduce two operators Si\ijkl) = \kjil), 

S2\ijkl) = \ilkj). They define three orthogonal subspaces of Hcomp ® T~Lcomp, 

Ua IIV') : 5i52|^) = -|^)} = span{|zjfcO - \kU])], (38) 
■H\ := (1^) : = ^alV-) = ^IV-)} = span{|iifcO + \klij) - \k31l) - \ilkj)], (39) 

nl := {H') : = S2W = IV)} = span{|zjfcO + \khj) + \kjil) + \ilkj)], (40) 



with V.A^B'HgiB'H'g — 'Hcomp®'Hcomp- Applying L to the vectors spanning the subspaces written explicitly in ([38l) - p0)) 
and using ([M)) and (|37l) . we find that ^-LA^ Hg and are eigenspaces of L with the eigenvalues, respectively, 

4 2 4 2 4 
Aa = 2 ^, Al = 2 ^, A^, = 2 + ^. (41) 

The largest eigenvalue Imax equals to A5 , the space Hg is thus in the kernel of the operator A and does not contribute 
Kraus operators to the generalized concurrence ca (|19p . The Kraus operators constructed from the vectors in Ha 
take the form T ~ \ij}{kl\ — \kl){ij\ [cf. (IT7)) ]. where \ij) = \i) (g) \j) etc. They are thus antisymmetric and, as explained 
in Section ITVl are also irrelevant to ca- We are left with the only ingredients given by the subspace Hg, i.e. the 
Kraus operators Tij^i = \ij){kl\ + \kl){ij\ — \kj){il\ — \il){kj\. Finally thus, in the definition we can take as A the 
projection on Hg, which is a subspace of Hcomp ® 'Hcomp = 'H®'H®'H®'H obtained by antisymmetrizing the first 
with the third factor as well as the second with the fourth ones. Calculations for T-Li 7^ 7^2 are only slightly more 
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complicated and give the same result, i.e. A as the projection on the subspace of Hi (E) H2 <E} Hi H2 obtained by 
separate antisymmetrizations of Hi factors (i.e. the first and third ones) and H2 (the second and fourth factor), which 
coincides with the results of [1, 01 . 

In the case of two identical particles with TV-dimensional single-particle spaces H the relevant representation is that 
of SU{N) on Hcomp — H AH or Hcomp = HV H. The corresponding representation of sl^iC) is given by 

Tr{A) ^ A (E) I + I (E A. (42) 

As for the case of distinguishable particles we start with the calculation of the Casimir operator C2. To shorten the 
considerations we may calculate C2 on H E) H for the representation (|42p . for which it is not proportional to the 
identity since the representation on the full tensor product is not irreducible. However, a further reduction to the 
irreducible invariant subspaces HAH and HV H will lead to the desired results. Again from ^ and (H^ we have 



C2 = ^ {XaX^a «)/ + /«) XaX^a + 2X„ E) X^a + X^^Xa + X^^X^ + 2X_„ ® X^) 

a>0 
N-1 

+ ^{Hf ®I + I®H'f + 2H, (g) H,). (43) 



Using the explicit form of Xij and given by ([5^ and ([55)) with the help of ([M]) we obtain 

(44) 

where the upper sign is for the symmetric subspace, H \/ H, and the lower one for the antisymmetric one, H AH. 
As it should, in each of these subspaces the Casimir operator is proportional to the identity, and in order to simplify 
calculations it is convenient to consider the operator 

L' = N{L~C2<E>I-I(E>C2) 

r 

= 2iV^ (7r(X„)®7r(X_„) + 7r(X_„)(g)7r(X„)) +2iV^7r(iJ,)(8)7r(i7,)- (45) 

Substituting ([32]) and ([33]) we obtain after straightforward calculations 

L'\ijkl) = (1 - 6,k)\kjil) + (1 - 5u)\ljkt) + (1 - Sjk)\ikjl) + (1 - Sji)\ilkj} 
4 

'n 



+ S^k + Su+Sjk+Sji ) \ijkl), (46) 



for \ijkl) eH(E)H(E)H(g)H. 

We are now in the position where we have to specify further calculations to Hcomp — Hy H (bosons) or He 
Hah (fermions). 



A. Bosons 



We choose a basis in Hcomp = H\/ H consisting of vectors 

\^^) = \i)E>\t), j = l,...,N (47) 
= + I,J = 1,...,N, i^j, (48) 

and split the space Hcom.p®Hcomp ~HyH®HyH into two parts, Hg which is symmetric with respect to interchange 
of the two copies of Hcomp and H^ which is antisymmetric with respect to this interchange. The space H^ can be 
decomposed into invariant spaces of the operator L', 

1. Hi spanned by vectors E) , 

2. H2 spanned by vectors <8) IV'j) + IV'j) ^ IV'i) and E) i 7^ j, 

3. "Hf spanned by vectors \4>i) (g) li/ij-,) -I- \4>ij) (g) j-^i), i ^ j, 
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4. "Hf spanned by vectors \^pi) ^ \ip.jk) + \i>]k) \^pi) and IV'y) ^ \ipik) + li'ik) ^ IV'y), i J ^ k ^ i, 

5. "Hf spanned by vectors (g) + IV'fci) IV'u), IV'ifc) li^ji) + li'ik), and iV'ii) «> IV'jfc) + IV'jfe) 
with all J, j, k, I different. 

Using (j46p we can diagonalize i' in each of the above subspaces. In particular we have: 

1. In the subspace 71 f 

4 



4- — ) |V'^)®|V'^)=A+|V'^)®|V'^)• 



2. In the subspace ^2 



A consequent diagonalization of the 2x2 matrix, 



4 

2 2 4 



reveals that 



where 



-'hTi2 



= (4-|)l*l,)=A+|*l,>, 



(49) 

(50) 
(51) 

(52) 



(53) 
(54) 

(55) 
(56) 



3. In the subspace Ti^ 



4. In the subspace 

+ 2{\^,j) (g) IV'jfe) + \^ik) (E) \'4'^J)) 

L'(IV'.y) (g li'tk) + li'tk) (g IV-'y)) = 4(|V'») (g IV'jfc) + IV'jfc) «) li/-*)) 

Hence, like in 2. above, 

i'l*?,.) =(4-^)l%) = A+|v|/3^.,), 

'ijk) ^ V ^ IV/ I ''ijfc/ 



where 



l^lk) 



L'\^t,k) =(-2-^)|*l,)=A_|*y, 

-2|l/^,) (g) IV'jfe) - 2|V'ife) (g \i>^) + g) IV'^fc) + |V'*fc) «) l^/'zj)- 



(57) 



(58) 
(59) 



(60) 
(61) 

(62) 
(63) 
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5. In the subspace 



+ 2(|V'ii) (8) IV'j/c) 4 



N 



+ |V'jfe)«> 



The relevant matrix to diagonahze reads as 



' N 



2 2 
2 

2 

^ AT 



which leads to 



L'^jki) 

' ijkl ) 
ijkl ) 



I* 



ijkl . 



ijkl/ 



= -2 



)Kki) = ^-Kki)^ 

)\nki)-^-Hki), 



with 



I* 



ijkl I 



I* 



ijkl I 



(64) 



(65) 



(66) 



(67) 



(68) 
(69) 
(70) 



(71) 
(72) 
(73) 



To finish the calculations in the bosonic case we have to consider the antisymmetric part ^HVHA'HV'Hof 



the space 'Kc 



It decomposes into invariant subspaces 



1. spanned by vectors 1*^^-) = (g) IV'i) - IV'j) 8) \ipt), i 7^ j, 

2. nf spanned by vectors \^fj) = (g) iV'y) - <E) IV'i), « 7^ j, 

3. ni spanned by vectors \^ljk) = Hi) <^ IV-'jfc) - IV'ifc) ® I-^j) and l^f^k) = IV'ii) ® IV'jfc) - Ife) (8 IV'ij), i J ^ k, 



4. -Hf spanned by vectors 1*1^-^;) = IV-jj) IV'fez) - IV-fci) IV-ij), 1*?^^;) 
l^ijki) = IV'ii) 8) IV'jfc) - IV-'jfc) 8) IV'ii) with ah i,j, k, I different. 

From (|46p we obtain straightforwardly 



\^tk) <8) IV'ji) - IV'ji) 8) and 



(74) 



for aU the subspaces, i.e. for |*) = \^f), and l^'lj'f). 

Furthermore, we easily see that if we do not assume that all i,j,k,l are different in the definitions of l^'^j^.f) and 
(adopting the notation \tpii) = 2\ii)), we can express all vectors spanning the subspaces , z = 1, ... ,9, in 



,T.l,2,3x 
ijkl I 



I* 

terms of {9]}^) and |#i (for example, IV'.) ® |^^,,) + IV'y) ® = IMu.) etc.). 
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To sum it up, the Hilbert space Hcomp "SJ T-Lcomp = 'HyH(E)HyH splits into three eigenspaces of L', 

H?-span{|vl/i^.,,)}, (75) 
[see Equations (|49l),(l53|,(l53,(l5B,(l60|,(|^ and 1^], 

= spanll^-i^.,;), 1*2^.,,), 1^-3^.,,)}, (76) 

[see Equation ([7^]. and 

n^=sv&n{\^lki)A^%ki)], (77) 

for aU i,j,k,l, not necessarily different - see the remark above, [see Equations (iMt . (|56| . (|6T|) . ([63t . (|69| . (ITO)) . (|72t . 
anddZl]. 

They correspond to three different eigenvalues of L', respectively, A+ = 4— A = — and A+ = —2 — The 
subspaces (|75II77[) are also eigenspaces of i = jjL' + 2 (l ~ jl^- + ■^) /, and the corresponding eigenvalues of L are 
= 2— + A^ — 2 — ^ + jj, and A^ = 2 — The largest one, Imax — ^+ corresponds to the kernel of 
the operator A = ImaxI — L, and vectors in the subspace give antisymmetric Kraus operators T^. Hence as A we 
may take the projection onH^ and construct the relevant Kraus operators from vectors in this space. 



B. Fermions 



Calculations for fermions closely follow the bosonic case, so we omit most of details. In the space Hcomp —HAH 
we choose a basis = |i) <E) \j) — \j) (E) j}. As in the bosonic case we decompose Hcomp ® Hcomp into the 

symmetric Hg = H AHV H AH, and antisymmetric, H^ = H AH AH AH, parts. The symmetric part can be further 
split into invariant subspaces of L' 

1. Hf spanned by vectors 1$,^^) — \(f>ij) (g) for which 

i'l*l,)= (2-^)|<i>l,)=AV|$i,.) (78) 

2. H2 spanned by vectors ~ \(f>ij) (g) \(t)ik} + \(f>ik} (Xi j 7^ k for which again 

i'l*?^.)^ (^2-l)|ci>5)=A;|$5) (79) 



3. H^ spanned by vectors \(j>ij) (g) 10^,) + 10^;) (g) \(j),k) «) + (g and g) \(j)jk) + |0jfc) <E) 10*;) 

with all k, I different. 

Here we diagonalize a set of equations 

L' {\(t),j) ® + \(t>kl) (S) \(t>^J)) ^ -J^ilM (E)\(t)^J) (S\(t)kl)) 

+ 2{\(j)^k) g) \(l>jl) + \(t)jl) g) |0»fc)) 

+ 2{\(j>u)(E)\(t>Jk) + \<l>3k)<E)\(f>^l)) (80) 

i'(|</>zfc) g) + g) |</>zfc)) = 2{\(j)ki) (S) \(t>tj) + \(t>tj) ® \<j)ki)) 

- ^{\(t)^k) <E) \<l)jl) + \(t>jl) (E) \(|)^k)) 

- 2{\(j)u)(E)\<|)Jk) + \<f>jk)<E)\(|)^l)) (81) 

L'{\(l)a)(E)\(t)jk) + \4'jk)(E)\(t)rl)) = 2{\(j)kl) (S) \(t)rj) + \(t)^J) \(j)kl)) 

- 2{\<P,k) <E) + \(f>jl) g) \<p^k)) 

- ^{\^^l)<»\M + \<pJk)<S>\(|)^l)). (82) 
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producing 

L'\'^liki) = (-4 - ^) |<&U> = A'_|$l,fe,), (83) 
L'\^lki) = (2 - ^) l^lki) = AVI$?,«), (84) 
L'l^U) = (2-^)l*?,feJ = AV|a>3.,,), (85) 

with 

l^ijki) = -l^u) ^ \M - o \<l>ij) + \(t>ik) o \(t>3i) 

+ |<?!>ii) |<^!>^fc) + \(t>^l) ® \(l>jk) + \(l>jk) ® |<^!>^^) (86) 

l*»jfci) = |'^!>^J)'»|'^!>fc^) + |0fe^)<^|0^J) + |'/>^^)<^|0Jfe) + |0J7c)®|0^^), (87) 
l*?,-fe() = l<^y) ® l-Z-w) + l<^fci) ® \4>io) + \<i>ik) ® \<t>3l) + \(t>3l) ® l<^ifc)- (88) 

The antisymmetric part 'H^ splits into invariant subspaces 

1. Uf spanned by vectors |f ^^.j.) = 10^) ® \(t)ik) - \(j)ik) ® j ^ k, 

2. ni spanned by vectors |l>|^-^;) = |(/>y) ® |(/>fc() - O |(/>ij), l^^i) = l^ifc) ® I'/'j') - ® l^ife)- and l^^-^^,) = 
\<Pii) ® \<Pjk) - \<Pjk) ® \<Pii) with all k, I different. 

In both subspaces L' acts as the multiplication by A' = — 

As in the bosonic case, we can express all vectors spanning subspaces Hf, i = 1,...,5, using only I^Jj^f) and 

l^ifkf) (here we have {(pa) = 0). 

Therefore L' has three different eigenvalues A'_,_, A', and A'_ corresponding to three eigenspaces, 

= span{|$|.,;),|$3.^;)}, (89) 



= span{|$i^.,;), m,^}, m,,}}, (90) 

and 

W^ = span{|4.,;)}, (91) 

for all i, j, k, I, not necessarily different. These are also the eigenspaces of i = j^L' + 2 (l — — ■^) /. The subspace 
Ti^ corresponding to the largest eigenvalue of L, Imax = ^+ = 2 — is in the kernel of A, whereas 'H^ corresponding 
to the middle eigenvalue A^ = 2 — -j^ — -^ofL gives rise to antisymmetric Kraus operators. The operator A can be 
thus chosen as the projection ouT-LZ- 
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